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ABSTRACT 

A procedure for the design of a recursive digital notch- 
filter with one coefficient is proposed. This type of filter 
can be implemented readily with CCD (charge -coupled device) . 
The effect of coefficient on notch-gain, notch-width, and 
passband-ripple is examined. The computer implementation is 
tested for second, third, and sixth order notch- filters and 
the optimum coefficient is found based upon limited computer 
wordlength . 
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I . INTRODUCTION 



Digital filters are finding wide applications in many 
diverse areas, such as speech processing, picture process- 
ing, sonar and radar system, digital control system, etc. 

A very useful type of such filters is the notch-filter 
which attenuates highly a particular frequency component 
in the input signal while leaving nearby frequencies 
relatively unattenuated. 

A linear digital filter can be realized physically in 
terms of the three basic operations of addition, multipli- 
cation and unit delay. The advantages gained by realizing 
a digital filter rather than its continuous counterpart 
lie principally in its stability, accuracy, reliability 
and flexibility [Refs. 7 and 20]. Furthermore, its 
frequency characteristics can be varied without altering 
the hardware simply by changing the stored set of filter 
coefficients. The disadvantages of digital approach lie 
in two extra design considerations. Firstly, the choice 
of wordlength of the computer and analog- to-digital con- 
verter. These problems do not arise in the corresponding 
continuous design. A design filter operates on an input- 
sampled signal by means of a computational algorithm. It 
can be simulated on a general purpose computer or can be 
constructed with special purpose digital hardware, such 
as CCD (charge-coupled device) implementation; this part 
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is beyond the the scope of research. Due to computer 
limited wordlength, the round-off error in computer program 
is produced. This can cause the filter characteristics to 
deviate considerably from the ideal. 

This research attempts to present a survey and 
discussion of the effect of the finite wordlength on the 
accuracy of digital notch- filter . The design procedure of 
digital notch-filter was to propose rational functions of 
a discrete variable which, xvhen transformed into the fre- 
quency domain, have desired filtering characteristics. The 
design was proceeded directly in terms of the variable, z, 
associated with discrete time intervals rather than 
through the intermediate variables associated with con- 
tinuous time. 

The digital notch-filter with one coefficient was 
investigated and tested by implementing on computer IBM 
360/70 installed at the Naval Postgraduate School. 
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II. THE CHARACTERISTIC TRANSFER FUNCTION 



OF DIGITAL NOTCH- FILTER 
A. TRANSFER FUNCTION 

As the first step in the design of the digital notch- 
filter, one often uses the bilinear transformation to 
approximate the analog filters in discrete time. This 
transformation process is perhaps the most utilized approach 
relating the continuous and discrete domains. The approach 
in this paper differs by starting directly in the discrete 
domain v^^ithout the bilinear transformation. 

The approach is to uniformly sample the signal u(t) every 
T seconds, and use the resulting sequence as the input to a 
recursive digital filter governed by the general difference 
equation 

v(K) = b^ u(K) + u(K-l) + ... + b^^ u(K-m) 

- a^ v(K-l) - a 2 v(K-2) - ... - a^^ v(K-m) (2-1) 

where u(k) and v(k) denote the values of the filter’s input 
and output signals, respectively, at the K-th instant. 

For minimum noise generating, cost, and simplicity, one 
can reduce equation (2-1) to the minimum number of coeffi- 
cients such that (2-1) is still in the recursive form. 

v(K) = u(K) + bj^ u(K-m) - a^^ v(K-m) (2-2) 
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The above simple recursive equation of the recursive 



digital filter is derived from (2-1) by letting 




and a, =a»=...=a n=0 

12 m-1 

Then its corresponding transfer function has the form 

1 + b 

2 * 

For a notch-filter one wishes the filter characterized by 
(2-3) to effectively transmit all sinusoidal inputs with the 
exception of the Wq radiant/sec. sinusoid which it should 
completely reject. Unfortunately, it is impossible to 
implement an ideal notch-filter 




Figure (2-1) Gain frequency behavior of 
an ideal notch filter. 
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B. GAIN FACTOR OF A DIGITAL FILTER 



The poles and zeros of the transfer function (2-3) are 
computed as follows 

z + b =0 
m 

hence z = ^/-b = e^ 

m m 

and z^ + a =0 

m 

n=0 , 1 , 2 , . . . 

m / m , — j (2n+l) — 

z = /-a = /a e-^ ^ 

m m 

The zeros and poles are respectively given by 

z. = (2-4) 

and 

— ^ / — i(2n"*'l) — ^''7c^ 

p . = /a e** ^ -^m (2-5) 

m 

i=l , . . . ,m . 

From (2-4) and (2-5) the zeros and poles are all complex 
conjugate for m even and for m odd, one of the zeros and one 
of the poles will be negative real. 

Let m=3, the zero's and pole's pattern is plotted in 
Figure (2-2) . 
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Figure ( 2 - 2 ) Pole-zero pattern. 

where 

is complex conjugate of 
T is sampling frequency, 

Wq is a notch- frequency which the one in which the first 
zero and pole occur. 
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For m-3 it is given by 



z 



1 





Pi 





After multiplying the numerator and the denominator of (2-3) 
by and factoring, the following zero-pole format is 
obtained 



* * 

(z-z,) (z-z, ) ... (z-z-p) (z-z-p ) 

H(z) = i ^ ^ (2-6a) 

(z-p^) (z-Pj^ ) ... (z-p2jj^) 2 ^ ^ 



where m is even and given by 



m = 21 

a = 1, 2, ..., N. 



If m is odd, called m = 2fi-+l, the format becomes 



H(z) 



(z-z^) (z-z^ ) ... (z-Z2j^) (z-Z2^ ) (^~^ 2 Z+l^ 

s s (2-6b) 

(z-p^) (z.p^ ) ... (z-p2jj^) (.z-p^^ ) (z-P 25 ,^i) 



where Z 2 ^^.-^ and ^ 2 i+\ negative real. 

The condition for the system being stable is: 



2 |H I < “ (2-7) 

o )=0 
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and the steady-state response of any stable system with 
transfer function (2-6) to the sinusoidal input, 
u(K) = Sin Ko)T, is given by 

y(K) = |H Sin (KuT + 9) 

where 

H (e^""'^) = H (z) I ^ ^ gjo)T 

and 

9 = angle {H (e^^"^)} . 

T CO T 

The term |H (e^ )| is called the system gain factor since 

it characterizes the sinusoidal gain response property of 
system (2-6). Hence, the steady state sinusoidal response 
is readily obtained by evaluating the system transfer func- 
tion at z=e^ , where oj is the radian frequency of the input 
sinusoid and T is the sampling period. One can determine 
the system's gain and phase angle factor in a straight- 
forward manner by using a geometrical approach. This approach 
gives the advantage of relationship between the locations 
of zeros and poles of the system transfer function and its 
gain-factor behavior. In the design of digital filter, it 
will enable us to design a system whose gain-factor approxi- 
mates any desired sinusoidal discrimination behavior. 
Substituting z=e^*^^ into (2-6b), 
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c?o src ore 




K-SCflLE:--:G.00t>0i UNITS INCH- 
V-SCRLt:---S. OOE-Oi UNITS INCH. 



Figure (2-5) Frequency magnitude plot for m=3. 
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(e^“* -Pj_) (e^“* -pj^') . . . (e^“' -p2j) (e^“‘ -p^j’') (e^“‘ 



(2 



where 



Similarly 



m rr — 3 (2n+l) — 

z-= /b e'^^ ^m 
1 m 

i = 1 , 2 , . . . , 11 
n = 0 , 1 , . . . , 1 % 
m = 2Jl + 1 



^ 2 Jl + l 



m 



/b^ . 



m 



Pi 



■"/S- 

m 



— 

m 



let 



2£ + l 



m 



m 



a 



gj«T 



m 



/b” e 
m 



j(2n+l) 



TT 

m 



= ej"''' - 



"•/i- ej( 2 n*l) 

m 



!L 

m 



_ ^jcoT . m r — 

“ 2 i*l ■ ® 



P 2 JI+I) 

8 ) 
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m\!r 




I 

( 



f 



3 



21+1 



jo)T 



= e 



m 



m 



Hence (2-8) becomes 



^ T> CIt o. 






a 



21+1 



^ 21 ^ ^^ 21 * 



3 



(2-9) 



2Z+1 



It is obviously the system gain factor at radian 

frequency w is seen to depend on the distance by which the 

system transfer function's zeros and poles are located from 

i T 

the critical point e^^ . Thus, if one wishes to completely 
reject the sinusoid sin kwT, it is necessary that the trans- 
fer function have the radian frequency w increases from 
zero to infinitive, one sees that the zero of transfer 
function (2-3) will fall at the first of the zeros of 
[(2n+l)— ] radians for n=0 and repeats for n=l,2,... as in 
Fig. (2-3). 




Figure (2-3) Gain frequency beliavior of 
digital notch- filter . 
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By considering pole-zero pattern in Fig. (2-2) and the 

equation (2-9) one wishes to reject completely at 
tt 

z = e^ the distance must be equal to zero. 

One knows that 

m A " - j (2n"*"l) — 
a . = z - /b e-' ^ ''m 

1 m 



hence 



a. = . 

1 m 

From the above equation, 

a. = 0 implies b =1 

and for the gain of (2-9) is constant for all frequencies 

TT 

except one at z=e^ . One should choose the coeffi- 

cient, a , as close to unity as possible: the closer a is 
to unity the closer the notch-filter approaches the ideal 
notch-filter as shown in Fig. (2-1). 

Looking back in Fig. (2-2), the parameter A = which 

controls the notch-bandwidth, when a close to one, 

A close to zero and so does the notch-bandwidth. Ideally, 
it is desired to make A as small as possible, but in prac- 
tice considerations one is not able to make A arbitrarily 
small (A will be discussed in section D of this chapter) . 

The above considerations can be made clear by analyzing 
the transfer function (2-3). 
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C. ANALYSIS OF MAGNITUDE RESPONSE 

It is discussed here the magnitude response as a 

function of frequency and the coefficients a and b . 

mm 

The equation (2-3) is rewritten as follows: 

+ b 

HCz) = -ji 21 ( 2 - 10 ) 

z + 

m 

Let z=e^“^ substitute into (2-10), it is given by 



H(e^“’^) = 



eJnuDT + b^ 



j'mwT ^ 



m 



( 2 - 11 ) 



i X 

Using e*^ = cosx + jsinx, one gets 



. T, (Cos mwT + b ) + j sin mwT 

H(e^“^) = ^ 

(Cos mo)T + a ) + i sin mwT 



( 2 - 12 ) 



• rp 

The squared magnitude of H(e^^ ) will be 



. _ _ (Cos mwT + b )^ + sin^ mwT 

l„ , 103T. I 2 ^ m^ 

|H(e-’ )| = 2 2 

(Cos mcoT + a ) + sin mu)T 

m 



or 

1 + b ^ + 2b Cos mwT 

(2-13) 

1 + a + 2a Cos mojT 
m 

Take derivative of (2-13) with respect to w and let the 
numerator equal to zero 



H(e^“h I 
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- 2mb„T sin muT [1 + + 2a^ Cos muT] 

m ‘■mm ■■ 

+ 2ma T sin mwT [1 + b ^ + 2b Cos mcoT] = 0 
m ‘•mm ■' 

or 

(-2mb T) (1+a^^) sin mcoT + 2ma T (1+b sin mmT = 0 
^ m m"^ m m-^ 

Hence 

sin mmT [a„(l + b„^) - b„(l - a^^)] = 0 (2-14) 

Suppose the term in the bracket is not equal to zero. 

The frequency at which the squared magnitude is maximum or 
minimum will satisfy the equation (2-14) 

n TT 

sin mmT = 0 implies coT = — 

^ m 

n= 0 , 1 , 2 , . . . 

At n=0, one has o)T=0, the magnitude has the maximum 

7T 

value and at n=l, mT=— , the magnitude has the minimum value 

and m=— is defined as a notch- frequency ; when n=2,3,.... 
the magnitude function just repeats as in Fig. (2-4). 

If one wishes to get the constant gain except at frequency 

TT 

w=njjjTp , n=l,3,5,..., as in Fig. (2-3), the coefficients a^^^ 
and b„ are selected such that this result is achieved. 

To do that in equation (2-14), one lets the factor in the 
bracket always equal zero, i.e.. 
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Figure (2-4) Gain frequency behavior with some 

values of a and b . 

m m 



a^ (1+b^^) 
m m 



b (1+a "•) 
m ^ m 



0 



or 



a + ab^-b - ba^ = 0 
m mm m mm 



= a„b^ (a„-b_) 
mm m m m m 



(2-15) 



Hence a b =1 

m m 

We know that the system is stable when all the roots 
of characteristic polynomial or transfer function of 
system are inside the unit circle of z-plane. 

By that condition, the coefficient a^^ must be 

a„ < 1 . 
m 

As stated in section B, the coefficient b has the value of 

’ m 

one for completely rejecting at the notch- frequency . One 
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sees that with bj^=l and a^^ < 1 , the equality (2-15) cannot 
be satisfied. But we can approximate (2-15) by letting the 
coefficient a^^ very close to one. 

One sees that both statements by analyzing the transfer 
function (2-3) and considering the pole- zero pattern in 
Fig. (2-2) are appropriate. 

The transfer function (2-3) now becomes: 



where the subscript m is deleted and its squared magnitude 
is given by: 



frequency w for m=3 in Fig. (2-5). The half power frequen- 
cies are that at which the gain is 3dB down or 




(2-16) 



lH(ej“'^) I 2 



1 + a^ + 2a Cos mwT 



2+2 Cos mtoT 



(2-16a) 



The plot of |H(e^'^"'^)| for particular value of a vs radian 



|H(e^“’')| - i 



or 



1 + a^ + 2a Cos mwT 



2+2 Cos mwT 



1 

2 



(2-17) 



Solve (2-17) for w 



2 

4 + 4 Cos mo)T = 1 + a + 2a Cos mwT 



(4-2a) Cos mcoT = a^ - 3 
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Hence 



Cos mojT = 2a (2-18) 

There are two solutions and 0 J 2 T which are symmetric 

to wT=^ , the notch- frequency of filter, as illustrated in 
Fig. (2-5) 




m 



Figure (2-6) 



where B is a notch-bandwidth. 

The notch-bandwidth B will be 



B 'T'/' ^ 1^ — 1^*^ — 

^ = T(a). - Wt ) = 2 [— - - Cos (-j ^) J 

27t ^2 1 '■mm 4 - 2a ■* 



or 



D 4tt r r- "1 - 3 

B = — [tt - Cos (-{ ^ 

m ^4 - 2a 



)] 



(2-19) 
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I£ one inserts the parameter A=l-a into equation (2-19), 
the bandwidth will be a function of A 



B = 



4tt 

m 



[tt - Cos 



-1 (1 - A )^ - 3 . 
4 - 2(1 - A)- 



or 



B = — [tt - 
m 



Cos 



-1 A - 2A - 2- 
2A + 2 



( 2 - 20 ) 



From (2-19) and (2-20), the limit of the bandwidth B 
when the coefficient, a, approaches one (or A approaches 

zero) is zero. This is valid in theory; in fact, there is 

* 

some limit for the coefficient, a, (also for A) in practical 
view (this phenomenon will be seen in Chapter V) . 



D. LOWER BOUND OF A DUE TO FINITE PRECISION ARITHMETIC 

Because of the limitations caused by finite wordlength 
in all current digital computers, a real number must be 
represented by using a finite number of bits in either the 
fixed-point or floating-point form. The error introduced 
by such a representation is discussed in this section. Only 
binary arithmetic will be discussed here. 

In a fixed-point representation, assume a number x has 
been normalized such that |x| £ 1 can be represented in twos 
complement form as 

oo 

X = b„ + E b-2'^ 

0 i=l " 



where 



b^=l or 0 and b^ is the sign bit. 



24 



I 



To approximate x by a t-bit register length, one of two 
processes is generally utilized: rounding or chopping. 

In rounding, a 1 or a 0 is first added to the t-th bit, 
(x)^_^, according to whether the (t+l)-th bit, is 

1 or 0. Then only the first t bits of the result are kept. 
In chopping, those bits beyond the most significant t bits 
are simply dropped. Since the error introduced by chopping 
is larger than that introduced by rounding, chopping arith- 
metic is not commonly used. Therefore, the discussion 
below is limited to the rounding arithmetic process. 

Let be the t-bit representation of the number x. 

The error which will be introduced is defined by 

e = (x)^ - X 

where |e| 1 2 ^ (2-21a) 

It is known that if two t-bit fixed-point numbers are 
added, their sum still contains t bits, provided that there 
is not overflow. Therefore, under this assumption of no 
overflo^^^, fixed-point addition causes no error. On the 
other hand, the product of two t-bit numbers may have more 
than t-bits. Thus, rounding is necessary if only t bits 
are to be kept for future processing. 

In the case of a floating-point number, x is represented 
as 



x 



f2® 
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where £, a fraction between y and 1, is called the mantissa 
and e is a binary integer called the exponent. Let x denote 
the t-bit mantissa floating-point approximation of the 
number x. The difference between x and fl(x)=(x)^ is called 
the round-off error and depends on the size of x. It is 
therefore the best measure relative to x. Then it is clear 
that 



(x)^ = x(l + e) (2-21b) 

where the relative error e is bounded by -2 ^ ^ e £ 2”^ 
Unlike the fixed-point case, both addition and multiplica- 
tion in floating-point can introduce round-off error. 

From the preceding section, it is apparent that a 
digital notch filter with notch at can be implemented 

by means of a filter whose transfer function has a zero 
located at e^^O . In order that the gain of such a filter 
in the frequency domain is near one for frequencies suffi- 
ciently close to Wq, it is necessary that there exist 
properly located poles which tend to cancel the effects 
of these zeros. This can be accomplished by selecting an 
appropriate A; from equation (2-20) it is obvious that the 
notch of the filter is directly dependent on the parameter 
A. Ideally, it is desirable that A be made as small as 
possible in order to make the notch width small. 

There are practical considerations, however, which when 
taken into account, will not enable the filter designer to 
make A arbitrarily small (see Chapter V) . In the design 
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of a notch-filter of (2-16), one wishes the gain factor 
at to be zero. For the critical point e^*^^ sufficiently 

. TT 

outside the circle of center = e and of radius A, 

Fig. (2-6), the gain factor for all practical purposes is 
one . 




Figure (2-6) A pair of first pole and zero. 

This follows since, for z=e^^^ outside this circle, by 
selecting A small enough, the ratio of the product of the 
distances from zeros and poles to the critical point e^*^ 
is effectively a constant. Since all digital computers 
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are limited by finite wordlengths, the desired coefficients 
will be truncated and the coefficient inaccuracies will 
result as stated at the beginning of the D of Chapter II, 

The characteristic difference equation of a digital 
filter (2-16) is given by 

v(K) = u(K) + u(K-m) - av(K-m) (2-22) 

When one wishes to implement (2-22) on a digital computer, 
the coefficient, a, will be stored with error due to trunca- 
tion as (a+e) , where e is the round-off error in (2-21a,b). 
This will rise to a shift in the actual pole locations of 
the transfer function. 

For system stability, the condition placed on (a+e) is: 

a + e < 1 

or 1-a > e 

hence A > e (2-23a) 

If one uses a t-bit wordlength, (2-23) will become: 

A > (2-23b) 

In the implementation of digital notch-filter (2-22) 
with t-bit wordlength, one must select the parameter 
A > 2 ^ to guarantee that the notch-width is sufficiently 
narrow and to avoid the possibility of unbounded system 
growth. It is noted that if one counts on guard bit, the 
wordlength will be (t+1) bits. Given notch-width B, one 
can calculate A from (2-20) and therefore t. 
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III. STATE SPACE REPRESENTATIONS AND REALIZATIONS 



The transfer function (2-16) of an m-th order digital 
filter is the relationship between one input and one output. 
Only state and output equations with single input and single 
output will be considered. In the time invariant discrete 
system, the transfer function (2-16) can be described by 
a state variable representation of the form 

x(K+l) = A x(K) + B u(K) (3-la) 

v(K) = C x(K) + 6 u(K) (3-lb) 



where 

A is an m X m matrix, 

B is an m X 1 matrix 

C is a 1 X m matrix 

6 is a 1 X 1 matrix or scalar. 

Taking the z-transform of (3-la) and solving for X(z) yields 

X(z) = [I - Az‘^]‘^ B U(z) (3-2) 

where I is an m-dimension identity matrix. Then, taking 
the z-transform of equation (3-lb) and substituting (3-2) 
for X(z) yields 

V(z) = {C [I - Az’^]'^ B + 6} U(z) 

= H(z) U(z) (3-3) 
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where 



(3-4) 



H(z) = C [I - Az'^]'^ B + 6 



H(z) is the transfer function of system (3-1). The 
discrete system is considered, as stated above, for only 
one input and one output. Hence, the matrix parameters 
A, B, C and 6 have the forms: 



m X m 
A = 



m X 1 
B = 



11 



ml 



1 

0 



Im 



m X m 



(3-5) 



1 X m 
C = 



[C, 



• • 

m^ 



6 = d 
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There are many ways to realize the digital filter (3-1). 
Some of them are the canonic realizations. The term 
"canonic” is used in the sense of minimizing the number of 
operations required. This minimization implies a minimum 
number of coefficients of the difference equations, which 
is the source of noise. 

If one wants to describe the system (3-1) by a matrix 
form as : 



S = 




(3-7) 



where S is called characteristic matrix system. 

The matrix parameters A, B, C and 6 are respectively 
given as follows: 



' 0 

I • 

(m-1) X (m-1) 

I 

1 • 

I 0 

1 

0 



(3-6) 



0 
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0 . . . 0 ] 




To describe a system in the matrix 
equations will have the form: 

CK+1) = - a x^ (K) + 

X 2 CK+1) = x^ (K) 

X (K+1) = X , (K) 
m ^ ^ m-1 ^ ^ 

V (K) = (^) (K) 

The corresponding realization of (3-7) 



form, its state 



u (K) 



(3-7) 



+ - u (K) 
has the form 
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u(K) 






1 

a 



a-1 




Figure (3-1) Non-canonic realization o£ m-th order 

digital notch- filter . 

The realization in Fig. (3-1) is not canonic; the purpose 
of this paper is to use only one coefficient. There are many 
disadvantages in the above realization. Among these are 
the production of more noise, higher cost and greater network 
complexity. Hence, the matrix form (3-7) is a poor way to 
describe the system. This leads to the consideration of 
a search technique for the canonic realization as stated 
in section A. From this concept, one can realize the digital 
notch-filter (2-6) as follows: 
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v(K) 



Figure (3-2) 


Canonic realization of digital 






notch-filter . 


the corresponding 


state equations of realization 


listed below 






^1 


(K) 


= -a Xjj, (K+1) + u(K) 


^1 


(K+1) 


= x^ (K) 


X2 

• 


(K+1) 


= (K+1) 


• 

• 

m 


(K+1) 


= X , (K+1) 

m- 1 ^ ^ 


V 


(K) 


= Xj^ (K+1) + x^(K) 



(3-8) 



One sees the realization in Fig. (3-2); with this 
technique there is only one coefficient and it provides 
the minimum number of the coefficients of the m-th order 
recursive digital filter one can achieve. Therefore, in 
noise analysis, if one consider only the source of noise 
produced by the multiplications of the coefficients, it 
is not necessary to assume the sources of noise are 
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uncorrelated, since there is only one source of noise in 
the realization. By this advantage, in the design of high 
order of digital notch-filter it is not necessary to imple- 
ment by cascading first or second order sections as is 
often done at present. 

The "dual" canonic realization of Fig. (3-2) can be 
derived. The so-called "transposed realization" has been 
shown by Jackson [Ref. 4] to be given in the form 




of digital notch- filter . 



The corresponding equations are given by: 
(K) = -a u(K) + u(K) 

(K+1) = x^ (K) 

X 2 (K+1) = x^ (K+1) 



X (K+1) = X , (K+1) 

m ^ ^ m-1 ^ ^ 

V(K) = Xj^ (K+1) + u(K) 



(3-9) 
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In the state equations (3-8) and (3-9), the value of 
the coefficient, a, is taken between 0 and 1. It is known 
that the notch-gain and notch-width of a digital notch- 
filter are functions of the coefficient a. The question 
which arises here is: what is the optimum value of the 

coefficient a which satisfies the desired characteristics 
of the notch- filters and how large is the ripple in their 
passband? 

These questions are investigated in the following 
chapters by simulating the digital notch-filter with the 
digital computer IBM 360/70. One can see the difference 
between theory and practice in Chapter III. 

For example, the gain in equation (2-9) at notch 
frequency, in theory, is a decreasing function when the 
coefficient a is increasing. However, in practice, some 
values of coefficient a make the notch-gain change non- 
monotonically because of round-off noise in iteration 
computations in the program. 
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IV. IMPLEMENTATION 



A. COMPUTER IMPLEMENTATION 

In this chapter, a digital notch-filter which is 
characterized by the state equations (3-8) is implemented 
on the IBM 360/70. Plots of the gain-frequency response 
and the time response are provided. 

The input signal is a sine-wave, i.e., a signal of 
single frequency, which sampled every T seconds. The 
corresponding gain for each input sine-wave is computed 
when the steady state is reached, by the ratio of the value 
of the peak input sine-wave divided by the value of peak 
output response. Then the input sine-wave is scanned over 
the frequency range of the notch-filter and the gain- 
frequency response is plotted. It is noted that the notch 
of the filter falls at the first zero of filter; that is, 

. V 

^m 

at Zj^ = e . If one desires the notch-filter at Wq 
radians, the sampling time T is computed as follows 



or 






7T 

m 



Hence 



T 



ir 

nuo 



0 
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or 



T 



(4-1) 



1 



0 

where ^ ^ notch- frequency in hertz. 

The following sections contain implementations of 
second, third, and sixth order digital notch- filters on 
the IBM 360 computer using single precision arithmetic. 

1 . Second Order Digital Notch-Filter 

The results of m-th order digital notch-filter 
in Chapters II and III are applied. The second order 
digital notch-filter is characterized by the transfer 
function 

1 + z‘2 

H(z) = ± ^ (4-2) 

1 + az 

and it is realized as follows: 




v(K) 



Figure (4-1) Canonic realization of second order 
digital notch-filter. 
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The corresponding state equations are given by: 



^1 


(K) = 


-a X 


2 (K^) - 


^1 


(K+1) = 


^1 


(K) 


^2 


(K+1) = 


^1 


(K+1) 


V 


(K) 


^2 


(K+1) + X 



u (K) 



1 (K) 



(4-3) 



I£ one wishes a notch frequency at £q = 60 hz, the 
equation (4-1) gives the sampling frequency T as follows 



T 



J. ^ 1 

2mfQ ~ 2 X 2 X 60 



4.166666667 x 10 



The figures from (4-2) to (4-10) illustrate the 
gain frequency responses in dB with different coefficients 
of one thousand iterations. It is noted that at one thousand 
iterations the steady state of frequency response is not 
yet reached. These plots merely depict the behavior of 
the gain- frequency response with changing coefficients. 

The result of the simulating program gives the steady 
state of frequency response at about six thousand iterations. 
They are shown in Figs. (4-2) to (4-17) corresponding to 
some typical coefficients. And the notch-gains, notch- 
widths, and maximum passband-ripple in dB of the steady 
state gain- frequency response are tabulated in Table (4-1). 
Fig. (4-11) shows the plot of notch-width vs. coefficient 
in theory by evaluating the equation (2-19). 
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The time response of a second order digital notch- 
filter when the input is a sinusoidal signal of frequency- 
60 Hz, i.e., at the notch frequency given below. In theory, 
the notch-filter with 60 Hz notch will immediately reject 
the above input signal, but in practice the time response 
sinusoidal signal will be damped in some finite time, small 
or large according to the coefficient and order of notch- 
filters one utilizes. Figures (4-18) to (4-21) show the 
time response of second order 60 Hz notch- filters with 60 Hz 
sinusoidal input signal applied. 



40 






I > 

•-f 



|H(z)| 




- 86.25 dB 



• Oi'c'-fOi UNITS INCH. 

.OOC-i-lII UNITS INCH. 

Figure (4-2) Frequency response of second order 
digital notch-filter. a=0. 90714359 
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Figure (4-3) Second order 

digital notch-filter. 
a=0 .969901345 
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Figure (4-4) Second order digital notch-filter. 



a=0.979 
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K-SCRLt:---2.00f + 0i UNITS INCH. 
Y-SCHLE--i.OOt>Oi UNITS INCH. 

Figure (4-5) Second order digital notch-filter. 

a=0.98 
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Figure (4-7) Second order digital notch-filter. 

V 

a=0.9815 
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X-SCfiLE::=2.00fc'40i, UNITS INCH. 
'4-SCRLE:=2.Q0E+Qi UNITS INCH. 



Figure (4-8) Second order digital notch-filter. 

a=0 . 982 
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Figure (4-9) Second order digital notch-filter. 

a=0. 98298407 
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Figuer (4-10) Second order digital notch- £i Iter . 

a=0. 99845153 
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Figure (4-11) Plot of notch-width vs coefficient, 

in theory, of equation (2-19). 
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K-SCPLf>-2.0Dr+0i UNITS INCH. 

Y-SCPLf>1.00f*0i UNITS INCH. 

Figure (4-13) Steady-state gain frequency response 
of second order notch-filter with 6000 
iterations. a=0.985 
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Figure (4-14) Steady-state frequency response 

of second order notch-filter. 
a=0.99 
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Figure (4-15) Steady-state gain frequency response' 

of second order notch-filter. 
a=0.995 
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Figure (4-16) Steady-state gain frequency response 

of second order notch-filter. 
a=0 .999 
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Figure (4-17) Steady- state gain frequency response 

of second order notch- filter . 
a=0.9999 
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K-SCRLE=4.00E+02 UNITS INCH. 
Y-SCRLE=^5.00E-01 UNITS INCH. 
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X-SCnLt->i.00t>02 UNITS INCH. 
V-SCRLP:=5.0GE--0i UNITS INCH. 





Jh 




<D 


in 




C 


?-» 


o 


O 


•H 








nj 


C 




o • 


0) 


U Jh 




0) <D 


•H 


c/) 


o4h 


tH 
M-l *H 


^ O 


O 4h 


• 


<D X 


o 


C/) (J 


c 


C 4-J O 

o o ^ 
Cl. C! oo 
cn 

0) 1— 1 CN 


lD 

CTJ 


?-• 03 cx) 


C3 






Time 
digi 
a=0 . 




Sic cco oto- 



S9 



K-SCRLE==1.00E+02 UNITS INCH. 
V-SCRLE----5.00E-01 UNITS INCH. 



N 



to 

C 

o 



5 h 

O 

4 -^ 



M-i 
r- O 

o 

o 



Jh 

<D 

Jh 

O 



O 

c 



nd 

C • 

O Jh 
O O 
<D 4 -i 
to rH 




•H 
M-i M-! 

O I 

o ^ 
to O to 
^ 4 ^ in 
O O 1 H 
Cl, {=; LO 

to ^ 

<D rH 00 

Jh 03 Cn 
4 -J 

d) *H • 
B bOO 
•H *H II 
H TJ c 3 



(NJ 






<D 

C 3 

DO 

•H 

tL| 






oiu ^cc ai: mi- cio- 



60 



«-SCnLE-i.00e+Q2 UNITS INCH. 

V-SCRLe:^S.Q0f-0i units inch. 
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It is noted that the passband ripple is large for 
frequencies very close to the notch- frequency . Further, 
about one hertz from the notch, all passband starts less 
than 0.09 dB in magnitude and smoother for going farther 
the notch. 

2 . Third Order Digital Notch-Filter 

The third order digital notch-filter is charac- 
terized by the transfer function 

1 + z'^ 

H(z) = ^ ^ V (4-4) 

1 + az 

and realized as follows 




Figure (4-22) Canonic realization of third-order 
digital notch- filter . 
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The corresponding state equations are derived as 



^1 


(K) = 


- a X 2 (K+1) + u (K) 


^1 


(K+1) = 


(K) 


^2 


(K+ 1 ) = 


(K+ 1 ) 


^3 


(K+1) = 


X 2 (K+1) 


V 


(K) 


X 3 (K+1) + x^ (K) 



and the sampling time required for 60 Hz-notch is given as 

T - 2^ = 2 X U 60 = 2-”77777S x lO'^ sec. 

Similar to the results of the second order case, 
the frequency responses of a third order digital notch- 
filter are shown in Figs. (4-23) to (4-31) with one thousand 
iterations. The steady state of third order is reached in 
about ten thousand iterations. The steady-state gain 
frequency responses are illustrated in Figs. (4-32) through 
(4-37) . 

The time responses of third order digital notch- 
filter are shown in Figs. (4-38) to (4-41). Similarly, 

Table (4-1) of second order. Table (4-2) of third order 
is tabulated. 
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Figure (4-23) Frequency 
response of third order 
notch-filter. 
a=0. 90714359 



70.4 dB 



X-5CFLe;----2.00f.40l 'UNITS INCH. 

v-scRLT:.-=i.oot;+Oi units inch. 



64 




X-SCRLt->2.00E+0i UNITS INCH. 

'^-SCRLf.-i.OOt'.+Ol UNITS INCH. 

Figure (4-24) Frequency response o£ third order 

notch- filtre . 
a=0. 96901345 
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Figure (4-25) Frequency response o£ third order 

notch-filter . 
a=0.979 
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Figure (4-26) Frequency response of third order 

notch-filter. 

a=0.98 
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Figure (4-27) Frequency response of third order 

notch- filter . 
a=0.981 
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Figure (4-28) Frequency response of third order 

notch-filter. 

a=0.9815 
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Figure (4-29) Frequency response of third order 

notch-filter. 

a=0.982 
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Figure (4-30) Frequency response o£ third order 

notch-filter. 
a=0. 98298407 
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Figure (4-31) Frequency response of third order 

notch-filter. 
a=0. 99845153 
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Figure (4-32) Steady-state frequency response of third 

order notch- filter . a=0.95 
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Figure (4-33) Steady-state frequency response of 

third order notch-filter. a=0.985 
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Figure (4-34) Steady-state frequency response of third 

notch-filter. a=0.99 
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Figure (4-35) Steady-state frequency response of 

third order notch- filter . 
a=0.995 
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Figure (4-36) Staedy-state frequency response of 

third order notch-filter. 
a=0.999 



7 & 




K-jC.RLF-- 2.Ci0e>0i UNITS INCH. 

y-dcrlf-'U.oof-oi units inch. 

Figure (4-37) Staedy-state frequency response 

of third order notch- filter . 
a=0 .9999 
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3 . Sixth Order Digital Notch-Filter 

The sixth order digital notch-filter is 
characterized by the transfer function 



H(z) 



1 + 
1 + 




(4-6) 



One often realizes the sixth order digital filter 
by cascading three sections of second order filters. This 
is done to reduce the sensitivity of pole positions due 
to truncation because of round-off error of iteration 
computations in the simulation program. 

The sixth order digital notch-filter characterized 
by transfer function (4-6) is discussed here. There is 
only one source of error; the noise produced in the 
computer implementation is small compared with the general 
form of sixth order digital notch-filter of thirteen 
coefficients. Therefore, one can implement the sixth 
order digital notch-filter (4-6) by direct canonic reali- 
zation as given by Fig. (4-42). 
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sixth order digital notch-filter. 



The corresponding state equations are given by: 



^1 


(K) = ■ 


‘ a X 


6 


+ u (K) 


^1 


(K+1) = 


^1 


(K) 




^2 


(K+1) = 


^1 


(K+1) 




^3 


(K+1) = 


^2 


(K+1) 




^4 


(K+1) = 


^3 


(K+1) 




^5 


(K+1) = 


^4 


(K+1) 




^6 


(K+1) = 


^5 


(K+1) 




v(K) 


^6 


/ — \ 
+ 
V / 

+ 


x^ (K) 



One still requires the notch- frequency at = 60 Hz; 
the equation (4-1) gives the sampling frequency T as follows 



rj, ^ ^ 1 

Zmf^ 2 X 6 X 60 



1.388888889 x lO’^ sec. 



It is noted that the sampling time for sixth order is one 
half of that for the third order case. 

Figures (4-43) to (4-46) illustrate the time response 
of a sixth order digital notch-filter with input sinewave of 
60 Hz frequency. 

Figures (4-47) to (4-49) are frequency response plots 
of one thousand, two thousand and two thousand and five 
hundred iterations, respectively, corresponding coefficient 
a = 0.98298407. 

It is noted that the sixth order digital notch-filter 
requires about fifteen thousand iterations to reach the 
steady-state. The steady-state gain frequency responses are 
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illustrated in Figs, (4-50) to (4-55) and some numerical 
values are tabulated in Table (4-4) . 
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Figure (4*47) Sixth order digital notch-filter 

at 60 Hz notch and 1000 iterations. 
a=0. 98298407 
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Figure (4-48) Sixth order digital notch-filter 

at 60 Hz notch and 2000 iterations. 
a=0. 98298407 
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Figure (4-49) Sixth order digital notch-filter 

at 60 Hz notch and 2500 iterations. 
a=0. 98298407 
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Figure (4-50) Steady-state 
frequency response of 
sixth order notch-filter. 
a=0.95 
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Figure (4-51) Steady- state 
frequency response of sixth 
order notch- filter . ’ 

a=0 .985 



. ) 
c* 




- 70.2 dB 






n 

i:us 

n 

U ly \ . 



^01 



n 



uNIiS INCH^ 
UNITS \\<: h > 



9j6 




I 









Hz 



- 67 dB 



Figure (4-52) 



Staedy-state frequency response 
of sixth order notch-filter. 
a=0 . 99 
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Figure (4-53) Steady-state frequency response 

of sixth order notch- filter . 
a=0.995 
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Figure (4-54) Staedy-state frequency response 

of sixth order notch-filter. 
a=0.999 







■ : 





I 





- 2.2 dB 



K-Sr.RL"'- 2. DOe' -^01. JNIT5 INCH. 
^•-3rRL"--2j.022' 0.1 UNITS 1N2.H. 



Figure (4-55) Steady-state frequency response 

of sixth order notch-filter. 
a=0.9999 
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One is lead to attempt to realize the sixth order 
digital notch-filter by cascading three sections of second 
order as follows: 




Figure (4-56) Realization of sixth order digital 
notch-filter by cascade. 



The steady-state gain frequency response is plotted 
in Fig. (4-57) and (4-58) for coefficient a=0.99 and a=0.995, 
respectively. 

Compare with the results of direct realization of 
sixth order notch-filter in Figs. (4-52) and (4-53); the 
notch-gain in realization by cascading is deeper than in 
direct realization. But in passband it does not improve 
anything. It seems worse than direct form by introducing 
more noise; one can see gain frequency response behavior 
for coefficient a=0. 9999999 between them in Figs. (4-59) 
and (4-60) . 

It proves that the result is contrary with the 
statement of many current lectures: the pole-position of 

direct form is sensible than cascade form in realizing 
high order digital filter [Ref. 1] . The result in 
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% 




Figs. (4-59) and (4-60) states that: the pole-position 

in cascade form is much more sensible than direct form in 
realizing sixth order digital notch-filter characterized 
by the transfer function (2-16) . 
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Figure (4-57) Frequency response of sixth order 

notch-filter by cascade. 
a=0.99 
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Figure (4-58) Frequency response of sixth order 

notch-filter by cascade. 
a=0.995 
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Figure (4-59) Frequency response of sixth order 

notch-filter by direct form. 
a=0. 9999999 
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Figure (4-60) Frequency response of sixth order 

notch-filter by cascade form. 
a=0. 9999999 
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4. Observation and Discussion 



In theory, by inspection, the magnitude response 
in equation (2-9) and the notch-width in equation (2-19) 
vs. coefficient of m-th order digital notch-filter charac- 
terized by transfer function (2-16) , one sees that when 
coefficient, a, approaches one, the magnitude of notch-gain 
and notch-width will be smaller [see Fig. (4-11)]. From 
the magnitude response plots, it is seen that theory and 
practice are in agreement whenever the steady-state output 
response is reached. This can be seen in the steady-state 
frequency responses in Figs. (4-12) to (4-17) of second 
order. Figs. (4-32) to (4-37) of third order and Figs. (4-50) 
to (4-55) of sixth order. Some of numerical results of 
steady-state response are tabulated in Tables (4-1), (4-2) 
and (4-3) . 

For the case of intermediate state (e.g., state of 
one thousand iterations- -the steady-state responses are not 
yet reached), the theory and practice appear to disagree. 

One can see this in frequency response plots in Figs. (4-2) 
to (4-10) of second order. These phenomena are a result 
of the random errors produced by truncating coefficients 
and round-off error in iteration computations in the simu- 
lating program due to finite register length of the computer. 
In this case, it is due to limitations of single precision 
arithmetic on IBM 360. A second explanation of the apparent 
deviation from the theoretical stems from error in computa- 
tion of notch-gain in computer program: 
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Gain,„ = 20 log,„ Cl ^^!; °l I) 

AB ^10 'Peak of output*'^ 



The peak of input is always one, but in the computer program, 
the peak of input and output are picked simultaneously and 
it is difficult to choose the sampling time such that it 
exactly corresponds to input peak; this is demonstrated in 
Fig. C4-56) 




Figure (4- 56) . 

where the sampling time T is fixed for prescribed magnitude 
response of desired notch- filter . 

Suppose : 

the peak of input 0.9 
the peak of output 0.0002 
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The corresponding notch-gain will be: 

GAIN = - 73.06 dB. 

Now let the peak close to one; suppose: 
the peak of input 0.99 
the peak of output 0.0001 

GAIN = - 79.91 dB. 

The error thus realized approximately 6 dB for 
the notch-gain. By the difficulty of picking exactly the 
peak of input as well as the peak of output, the actual 
notch-gains are certainly somewhat deeper than the values 
listed in Tables (4-1), (4-2), and (4-3). Also, the actual 
passband ripples should be smoother. 

It is known that the second order needs around 
six thousand iterations, the third order needs around ten 
thousand iterations and the sixth order needs around fifteen 
thousand iterations to the steady- state . It is seen by 
considering the time-respond plots of second, third and 
sixth order. It is obviously true since higher order in 
direct realization requires larger number of delays. 

It is noted that the passband-ripple is large 
around notch- frequency and smoother when going further. 

The results show that the higher order of digital 
notch-filter improves the passband-ripple (i.e., reduced 
ripple) and notch-gain [i.e., deeper notch. See Tables 
(4-2), (4-2), and (4-3)] and steady-state frequency response 
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plots in Figs. (4-50) to (4-55) of second, third and 
sixth order, respectively. 

The sixth order notch-filter requires about four 
minutes to plot one steady-state gain frequency response 
using fifteen thousand iterations. 

One will see that the results in this chapter will 
give the easy way to work in optimum problems due to finite 
register length on computer in the next chapter. 

B. CCD IMPLEMENTATION* 

The charge-coupled device (CCD) was invented at the 
Bell Telephone Laboratory by Drs . Boyle and Smith in 1970. 
This new development has been receiving much attention 
by industry because of its important properties- -low power 
consumption, low noise, and delay capability. Many success- 
ful applications have been investigated such as digital 
signal processing especially in design of recursive digital 
filter. More details of CCD can be seen in many current 
lectures. The CCD hardware implementation of recursive 
digital notch-filter is beyond this research. The following 
is the sample of block diagram of second order CCD digital 
notch- filter : 



*”Charge-coupled Device for Analog Signal Processing-- 
Recursive Filters Study,” Engineer's Thesis by Varachai 
Imsa-ad, Naval Postgraduate School, December 1974. 



Ill 



u(K) 



v(K) 




Figure (4-61) Block diagram of a second order 
CCD digital notch- filter . 
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V. OPTIMUM COEFFICIENT 



In this chapter, a practical method is described for 
choosing the coefficients of a digital filter to meet 
arbitrary specifications of the magnitude characteristic. 

As shown above, all special purpose computers have 
finite wordlengths which is a major source of noise in the 
implementation of digital filters. This is because of the 
necessity to truncation coefficient length and also because 
of rounding in the arithmetic operations, especially in 
the iterative multiplications. This kind of noise deter- 
mines the signal- to-noise ratio of the system. In some 
cases, the limitation of wordlength of coefficients may 
lead to a system which does not satisfy the original 
conditions . 

Especially, in the design of recursive digital notch- 
filters, one has to select the coefficients such that they 
satisfy the prescribed characteristic of notch-gain, notch- 
width, and the smallest passband ripple. 

The m-th order digital notch-filter is characterized 
by the transfer function (2-16) . It is seen that the notch- 
gain, the notch-width, and the passband ripple are functions 
of a single parameter- - the coefficient "a" between 0 and 1. 
It is possible to construct a performance function of 
single parameter in the optimization problem. 
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A. SEARCH PROCEDURE 



In the design of an m-th order digital notch- filter , 
if one wishes the prescribed magnitude having the follow- 
ing characteristics; 

(i) Notch-gain is at least - 60 dB. 

(ii) Notch-width is at least one hertz, 
any coefficient, a, in )0,1( satisfies conditions (i) and 
(ii) , and that value which produces the smoothest passband 
is the optimum value. 

In the method of searching, first, one selects the 
subinterval of )0,1( such that any coefficient belonging 
to this subinterval will satisfy the conditions (i) and 
(ii) ; second, searching is accomplished for the minimum 
error function of single variable for one of the 
subintervals . 

The subinterval can be found by computer experiment. 
Suppose for sixth order notch- filter one can consider 
Table (4-3); call (a^ , a 2 ) this subinterval: 

ai a^ 

/ t€ S/ ,<///>'/ 

As stated above, any coefficient, a, belonging to 
interval (a^^ , a 2 ) will satisfy the conditions (i) and (ii) . 

For sixth order, the values of a.^ and a 2 are given as 
[see Table (4- 3) ] : 

a^ = 0.98 and a 2 = 0.995 . 
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For third order [see Table (4*2]: 

= 0.98 and a^ = 0.99 . 

It is noted that the ideal (errorless) magnitude 
frequency response is equal to zero dB over the passband 
of the notch- filter . Due to the limited wordlength of 
computers, the actual magnitude frequency response will 
deviate from zero dB over the passband. If one defines: 

e (a, = 1 - H (a, w^) (5-1) 

as the error function corresponding to coefficient, a, and 
radian frequency of sinusoidal input signal. The error 
function e (a, w^) is defined in discrete domain which is 
considered as approximates the mean-squared error in con- 
tinuous domain of interval between 0 and 

“i 

e Ca, M.) - ^ S |H„ (a, - H (a, dco, CS-2) 

where H (a, u.) is ideal response and has magnitude equal 

00 

to one . 

Rewrite (5-2) as 

to . 

e (a, ( 0 .) = ^ ^ |1 - H (a, co. ) | ^ dco. (5-3) 

1 Wi 0 

the system is stable, e (a, w^) has to be convergent. 
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Let scan over the frequency range of the passband of 
notch- filter , the error performance or the total error 
function, is defined as 

^ 2 

E (a) = E e ^ (a , w . ) 
i=l ^ 



or 



E(a) = Z [1 - H Ca, w-)] 
i=l ^ 



(5-4) 



where 3.2_ — ^ — ^2 N is a finite integer number. 

The optimization problem turns out to be the least 
square method of single parameter. The optimum coefficient, 
a, is one that makes the error performance E(a) minimum. 

The following algorithms are described to search for the 
minimum of E(a). First, the interval (a^, a^) is subdivided 
into N subintervals such that 



aj^ e [a^ + 



(K-1) 



N 



K(a -a ) 

> =*1 ^ 



where K = 1 , 2 , . . . , 
Then compare the 
to successive value 



a 



K+1 



e 




N. 



values of 
of aj^ and a 

K(a2-a^) 

N ’ 



E (aj^) 

K+1’ 




and E(aj^^^^ 
where 



corresponding 



(K+l)(a2 

iT 



) 

-] 



If E(aj^) ^ E(aj^_|_2^) 

the next step will compare E Caj^+^) ^ ^^K+2^ store 

coefficient 
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Or if E (aj^) £ E , the next step will compare 

E 3.nd E store coefficient a^,; the execution 

continues until the last subinterval of (a^, a 2 ^ is reached, 
The final result will be printed out as the optimum 
coefficient corresponding to the smallest value of error 
performance E (a) for a^^ £ a ± 

The greater the value of N, the higher accuracy one 
can get but the price is paid by spending more computation 
time . 

The plots of error performance are illustrated in Fig. 
(5-1) of second, third, and sixth order digital notch- 
filters. One sees that the error performance is monotoni- 
cally decreasing with increasing coefficient, i.e., the 
greater the coefficient, the smoother the response in 
passband. This result can be extended for many order of 
digital notch-filter characterized by the transfer function 
(2-16). The optimum coefficient now is so easy to pick and 
it is the coefficient ^ 2 - 

for sixth order a 2 - 0.995 

for third order a 2 - 0.99. 



B. DISCUSSION AND CONCLUSION 



One sees that the error performance of single parameter 
is unimodal. The optimum problem of this kind of function 
has been accomplished by many methods well-known such as 
Fibonacci Search, Golden Section, etc. The Fibonacci Search 
is known as the most efficient in optimization technique of 
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K-SCRLF-2. OOf-01 UNITS INCH. 
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Figure (5-1) Plot of total deviation of second, 
third and sixth order notch-filter 
between a=0.5 and a=0.999. 

The region in circle is exaggerated 
in Fig. (5-2) . 
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unimodal objective function of single parameter. The method 
used here- -referred to by some authors as the uniform or 
factorial search- -compares formatly with the Fibonacci 
Search in that it takes more computation time than, but in 
this case of error-performance E(aj, in the author's 
opinion, is the most efficient one. In spite of the cost 
of computation time, the method is simple and the result 
is precise. 

The results carried out in this chapter agree with 
computer experiments in Chapter IV by considering the 
numerical value in Tables (4-1), (4-2), and (4-3). Higher 
orders of notch-filter caused deeper notch-gains and smoother 
response in passband. It is noted that the smoother pass- 
band is correct for coefficient close to one. It is seen 
in Fig. (5-2) for noise-performance in the region close 
to one. 

The value of the optimum coefficient, a, increases with 
the order of filter provided the conditions (i) and (ii) 
are fixed. The higher order of notch- filter , the closer a 
approaches one, the constraint is placed on a as mentioned 
in section D of Chapter II, due to the finite wordlength, 
e.g., the single precision on IBM 360/70 has 32-bit 
wordlength. The inequality (2-23) gives 

- 31 

1 - a2 = A > 2 
or 

a2 < 1 - 2"^^ ^ 0.999999999 (5.5) 
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Figure (5-2) Noise performance in the region 
of coefficient close to one. 
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The value (1-2 is the upper limit of coefficient, a, 

in theory. In practical view, the upper limit is smaller 
than (1-2 ) by the reason of error introduced in 

arithmetic operation and error in representation as stated 
above. Besides, in the implementation of notch- filter , 
the optimum coefficient corresponding the larger wordlength 
in register of computer will be closer to unity than 
the shorter one. For example, in implementation of sixth 
order of single precision on IBM 360, the upper limit of 
coefficient is 0.999, for coefficient exceeds this limit, 
the system begins to degenerate. It is the same limit for 
second and third order. One can see this phenomenon 
numerically in Tables (4-1), (4-2), and (4-3) and the 
behavior of changing in steady-state frequency response 
plots of second, third and sixth order in previous chapter. 

One sees that for a = 0.9999, the notch-gain less 
than 3 dB in magnitude and the passband ripple is unacceptable 

It is foolish to accept them as the notch- filters . 

31 

Furthermore, for a > (1 - 2 ) the system may become 

the oscillator due to error in rounding since it could 
make the poles of their transfer functions exactly on the 
unit circle of z-plane. 
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VI. CONCLUSION 



The filter design was used to design the digital notch- 
filters. The effect of order of filter on notch-gain, 
notch-width, and passband ripple was investigated. 

From results of this study, the simplest configuration 
of digital notch- filters was found; hence implementation 
will be economical in hardware. The notch-frequency, the 
notch-gain and the notch-width are very easily controllable 
by varying only one coefficient in realization as well as 
the agreement between the predicted and practice sampling 
time. This property is usable in an array of notch- filters 
to eliminate interference signals of variety of frequencies. 

It was seen that with selected coefficient the pass- 
band is very smooth (for sixth order less than 0.05 dB) 
and the notch-width is very narrow (less than one hertz) . 

The improvement can be achieved by increasing the order of 
filters provided the coefficient close enough to the optimum 
coefficient which was found by Search procedure. The opti- 
mum computer program was used and the result was precise. 
This program required large computer time. That is the 
reason some authors do not like it. 
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